Global Persistence Exponent for Critical Dynamics 
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A 'persistence exponent' is defined for nonequilibrium critical phenomena. It describes the 
probability, p(t) ~ t~ e , that the global order parameter has not changed sign in the time interval 
t following a quench to the critical point from a disordered state. This exponent is calculated in 
mean- field theory, in the n = oo limit of the O(n) model, to first order in e = 4 — d, and for the 
1-d Ising model. Numerical results are obtained for the 2-d Ising model. We argue that 9 is a new 
independent exponent. 



For many years it was believed that critical phenom- 
ena were characterized by a set of three critical ex- 
ponents, comprising two independent static exponents 
(other static exponents being related to these by scal- 
ing laws) and the dynamical exponent z. Then, quite 
recently, it was discovered that there is another dynam- 
ical exponent, the 'non-equilibrium' (or 'short-time') ex- 
ponent A, needed to describe two-time correlations in a 
system relaxing to the critical state from a disordered 
initial condition . It is natural to ask 'Are there any 
more independent critical exponents?'. In this Letter we 
propose such an exponent - the 'persistence exponent' 
9 associated with the probability, p(t) ~ t~ e , that the 
global order parameter has not changed sign in time t 
following a quench to the critical point. We calculate 
9 in mean-field theory, in the n = oo limit of the 0(n) 
model, to first order in e = 4— d (d = dimension of space) 
and for the d — 1 Ising model. In fact, it turns out that 
all these results satisfy the scaling law 9z = A— d+l — r)/2, 
which can be derived on the assumption that the dynam- 
ics of the global order parameter is a Markov process. 
We shall argue, however, that this process is in general 
non-Mar kovian, so that 9 is in general a new, non-trivial 
critical exponent. 

The persistence exponent 9 was first introduced in the 
context of the non-equilibrium coarsening dynamics of 
systems at zero temperature [[§[|. In that context it de- 
scribes the power-law decay, p(t) ~ t~ e , of the probabil- 
ity that the local order parameter </>(x) has not changed 
sign during the time interval t after the quench to T = 0. 
Equivalently, it gives the fraction of space in which the 
order parameter has not changed sign up to time t. More 
generally, one can consider the probability po (ti , £2) of no 
sign changes between t\ and t 2 . Scaling considerations 
suggest po(h,t 2 ) = f{h/t 2 ) ~ (h/t 2 ) e for t 2 > t x . 

Exact solutions for one-dimensional systems in- 
dicate that, in general, 9 is a new non-trivial exponent 
for coarsening dynamics. Recently, we have shown that 
even the diffusion equation exhibits a nontrivial persis- 
tence exponent, and have developed a rather accurate ap- 
proximate theory for this case H . The diffusion equation 
is itself a model of ordering dynamics, via the approxi- 



mate theory of Ohta, Jasnow and Kawasaki (OJK) [Q, 
and also describes, in its essential features, the ordering 
kinetics of the nonconserved 0(n) model in the large- n 
limit ||: The exponents 9 for these systems (OJK and 
large-n) are just those of the diffusion equation. 

In this Letter we introduce and calculate the anal- 
ogous exponent 9 for non-equilibrium critical dynam- 
ics. In this case however, one needs to consider the 
global, rather than the local order parameter. This is 
because individual degrees of freedom ('spins', say) are 
rapidly flipping so that the probability of not flipping 
in an interval t has an exponential tail. We shall see, 
however, that the probability for the global order pa- 
rameter not to have flipped indeed decays as a power 
law. One simplifying property of the global order pa- 
rameter is that, in the thermodynamic limit, it remains 
Gaussian at all finite times. This follows from the cen- 
tral limit theorem, on noting that the order-parameter 
field (/>(x, t) has a finite correlation length, L(t) ~ i 1 / 2 . 
If the system has a volume V 3> L(t) d , the appro- 
priate Gaussian variable is the k = Fourier compo- 
nent, <f>o(t) — [j d d x 0(x, t)]/y/V. From standard scaling, 
(4>l(t)) - Lit) 2 -*]. This follows from the k -> limit of 
the scaling form <0 k (t)0_ k (t)) = k-^ 2 -^G[kL(t)]. 

Our explicit results are derived from the Langevin 
equation for the vector order parameter cj> = (</>i, . . . , 4> n ): 



V 2 ^ - rfa - (u/n)4>i $ + C 
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(1) 



where £(x, t) is a Gaussian white noise with mean zero 
and correlator (f f (x, t)^ (x', t')) = 25 a £ d (x-x') 5(t-t'). 
[For a vector order parameter, we are defining pit) as the 
probability that a given component of the global order 
parameter, $(£) = J d d x <j>ix.,t), has not changed sign]. 

In mean-field theory, valid for d > 4, we set r = = u. 
Then the k = Fourier component (/>i(0,t) = $i(t)/V~V 
(where V is the volume) obeys the simple equation (sup- 
pressing the index i and the arguments) 



(2) 



indicating that <f> executes a simple random walk. The 
non-flipping probability p(t) is therefore just the proba- 
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bility that the random walker has not crossed the ori- 
gin up to time t. It is given by |j] p(t) ~ \(j>o\/\/i 
for large t, where 0o is the initial value of <p. Finally 
one has to average over |0o|- For a disordered initial 
condition, $(0) ~ y/V by the central limit theorem, so 
<f>o = $(0)/vl7 is Oil), the desired average is also O(l), 
and p(t) ~ 1/ We conclude that 8 — 1/2 in mean-field 
theory. 

Next we consider the large-n limit. Equation |l]) then 
simplifies to the self-consistent linear equation 

d t< /> = V 2 0-(r + U <0 2 »0 + £, (3) 

for each component. Defining a(t) = —r — u(0 2 ) and 
b(t) = £ a(t') dt 1 , (g) has the Fourier -space solution 

0(k, t) = 0(0, t) exp[b(t) - k 2 t] 

+ f dt'i(k,t')exp[b(t)-b(t')-k 2 (t-t')} . (4) 
Jo 

It is easy to show that the second term, involving the 
noise, dominates the first at large t |Q. Retaining only 
the second, computing (0 2 ), and defining g = exp(— 2b) 
leads to the equation 

d t g = 2rg + 4u [ dt' g(t') V cxp[-2fc 2 (t - t')} , (5) 
Jo k 

which can be solved by Laplace transformation. Setting 
r equal to its critical value, r c = — fc -2 gives 

g(s) = [s + 4u{J(0)- J(s)}}- 1 , (6) 
J( S )=Y,( s + 2k 2 )- 1 , (7) 

k 

from which one deduces g ~ s ( 2 ~ d )/ 2 for s — > 0, for 
2 < d < 4. Inverting the Laplace transform gives (with 
e = 4 - d) g(t) ~ <~ £//2 for t — > oo, whence b ~ (e/4) lni, 
and a(t) — e/4t. 

The large-n equation of motion (g) therefore reduces 

to 

d t 4> = {t/u)~4> + i (8) 

for the k = Fourier component of 4>- The final step is to 
eliminate the first term on the right by setting (f> = t^^ip, 
to give dtip = i _e / 4 £(i). Introducing the new time vari- 
able t = t x , this equation reduces to the random walk 
equation d T ip — r)(r), with r\ a Gaussian white noise, if 
one chooses x = (2 — e)/2. The final result is therefore 

P (t) - T-V2 = t( 2 - d )/ 4 , giving 

9={d-2)/A, 2<d<4(n = oo). (9) 

For d > 4, 9 sticks at the mean-field value of 1/2. 

Finally, we calculate to first order in e = 4 — d. This 
is most simply accomplished using the method of Wil- 
son ITqI. To order e the calculation can be carried out 



in d = 4, by expanding p(t) to first order in u, setting 
u equal to its Renormalization-Group (RG) fixed-point 
value, and exponentiating logarithms. 

The perturbative calculation of p(t) is in principle 
quite a difficult task. A systematic technique for per- 
forming the perturbation expansion was recently devel- 
oped by two of us |l^] in the general context of first- 
passage-time problems for non-Markov Gaussian pro- 
cesses. It amounts to expanding around the random walk 
(2) within a path-integral formulation of the problem. 
Since the global order-parameter <fr(t) remains Gaussian 
at all times (in the thermodynamic limit), this method 
is applicable. In the present work, however, we restrict 
ourselves to first order in e, for which the result can be 
obtained by elementary methods. The reason is that the 
dynamics of Q(t) remain Markov to this order, as we shall 
see. 

First we replace u/n in ([!]) by u, to conform to the 
conventional notation for e-expansions. To first order in 
u, one can as usual replace the nonlinear term utfii <f> 2 
in (|l|) by the linearized form (n+2)u(<f> 2 )(j)i. The required 
expression for (</> 2 ) can be evaluated at u — and r = 0, 
since r c is also 0(u). For this part of the calculation, 
therefore, we can use the mean-field result, 

Mt) = 4> k (0)eM-k 2 t)+ f dt'i k (t') exp[-fc 2 (t-t')], 

Jo 

(10) 

to give 

($) = A^exp(-2fc 2 <) JL (1 - exp(-2fc 2 t)) , 

k k 

(11) 

where we have used (</>k(O)0-k(O)) = A, appropriate to 
short-range spatial correlations in the initial state. 

To order u the critical point is r c = — (n + 2)it^ k k~ 2 . 
The effective Langevin equation for the k — mode, cor- 
rect to 0(u), is therefore (suppressing the index i and 
the momentum subscript) 

dt4> = (n + 2)u (Jp - a) exp(-2k 2 t) + £. (12) 

The k-integrals are now performed in d = 4. It is clear 
that the term involving the initial-condition correlator A 
is smaller (by one power of t) than that coming from the 
thermal noise, and may therefore be dropped, giving 

dt $=(n + 2)^4> + i (13) 

where K4 = I/871" 2 is the usual geometrical factor. Set- 
ting u equal to its RG fixed-point value |l(J u* — e/[(n + 
8)Ki] gives an equation identical to the large-n equation 
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(g), but with the replacement e -> [(n+2) /(n+8)]e. Mak- 
ing the same replacement in (^) , we deduce immediately 
that the exponent 9 is given by 



e + 0(e 2 ), 



(14) 



which agrees with ( p|) f or n — ► oo. For the Ising univer- 
sality class (n = 1), lfl|) becomes 6* = l/2-e/12 + 0(e 2 ). 

The final soluble limit we consider is the d = 1 Ising 
model with Glauber dynamics. For this model, the crit- 
ical point is at T = 0, so the 'critical' and 'strong cou- 
pling' fixed points coincide. The persistence probability 
p(t) for a single spin has been considered earlier Q. Very 
recently, it has been shown to decay as i~ 3 / 8 , with non- 
trivial results for general g-state Potts models j5j. The 
calculation of p(t) for the global magnetization M(t) is 
much simpler. At T — the dynamics is equivalent to 
a set of annihilating random walkers (the domain walls) . 
At each time step, each random walker moves indepen- 
dently one step to the left or right p2[ . Let the number 
of spins be N . Then the number of surviving walkers at 
time t is of order Nf- 1 / 2 H,0. The change in M(t) 
in one time step is therefore of order ^/Nt^ 1 ^ 4 , since the 
contributions from the walkers add incoherently. The 
k = Fourier component cf> = M/ y/~N therefore satisfies 
the Langevin equation (up to constants) 



dt4> = *" 1/4 C(<) 



(15) 



where is a Gaussian white noise, (£(£)£(£')) = C8(t — 
t') 1 and C is a constant. 

This can be reduced to the standard random-walk dy- 
namics through the change of variable t = r 2 . Equa- 
tion (|l5|) then reads dcft/dr = 2t 1 / 2 £(t 2 ) = r](r), where 



t](t) has correlator {f](T)r)(r')) = ACtS(t 
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2C5{t — r'), i.e. rj(t) is a Gaussian white noise. We con- 
clude that p(t) cx t- 1 ' 2 = t" 1 / 4 , i.e. 9 = 1/4 for this 
model. It is remarkable, but certainly coincidental, that 
the O(e) result gives this result exactly, on setting e = 3. 

At this point we note a simplifying feature of all the 
results presented so far, namely the underlying dynamics 
is a Gaussian Markov process in every case. This should 
be apparent from equations (||), (||), ([l3|), and (|l5|). For 
such cases one can derive (see below) a scaling law relat- 
ing 9 to other exponents, namely 



9z = \-d+ 1 — 77/2 



(16) 



where A describes the asymptotics of the two-time cor- 
relation function of the global order parameter at T c : 
($(h)4>(t 2 )) = tf- v)/z F(t 2 /h), with F(x) ~ x^/z 
for x — > 00. Using the known results rj = 0, z = 2, 
A = (3d - 4)/2 for n = 00 f§, 17 = 0, 2 = 2, A = 
d - [(n + 2)/(n + 8)]e/2 to 0(e) Jj], and 7? = 1, z = 2, 
A = 1 for the d = 1 Ising model |13| |, one can check that 
all of the results derived above satisfy this scaling law. 



Does this scaling law hold generally? We do not think 
so: we believe that $(i) is not a Markov process in 
general (though it is Gaussian), for the following rea- 
son. Consider the autocorrelation function for the k = 
mode, {4>{tx)4>{t 2 )) . We have seen that it has the scal- 
ing form rf?~ v)/ *F(t2/ti), with F(x) - x^ d ^/ z for x -> 

00 . Now construct the normalized autocorrelation func- 
tion a{hM) = mS^/mi) 2 ) 112 ^) 2 ) 1 ' 2 - This 
has the scaling form a{ti,t 2 ) = f{t\/t2), with f(x) ~ 
x {\-d+i-7]/2)/z f or x _j QO _ we introduce the new time 

variable T = hit, this becomes A(T U T 2 ) = g{T x - T 2 ), 

1. e. the process is a Gaussian stationary process in this 
time variable. Also the function g(T) has the asymptotic 
form g(T) - exp(-A|T|), with A = (A-d+l - r)/2)/z. 

Now if the process is Markovian, g(T) necessarily has 
this exponential form for all T, not just asymptotically 



large T |15|. Futhcrmore, the first-passage exponent 9 is 
then equal to A |) 15 17| ], which is the origin of the scaling 
law ( |16|) for Markov processes. Note that, in the original 
time variables, requiring g(T) to be a simple exponen- 
tial is equivalent to requiring that the scaling function 
fiti/t?) be a simple power of t\jt 2 for all t% > t±, not 
just t 2 > t\. 

So the question of whether $(t) is a Markov process 
comes down to the question of whether the scaling func- 
tion f(ti/t 2 ) [of the normalized autocorrelation function 
of 4>(t)] is a simple power law for all t 2 >t\. The known 
results for n = 00, 0(e), and the d = 1 Ising model 
satisfy this requirement. For the last of these, the re- 
duction to a random walk makes this transparent. In 
the other two cases, it is consequence of the simplicity of 
the one-loop nature of the calculations, which give simple 
powers. To 0(e 2 ), however, the structure of the 'water- 
melon' two-loop graph leads to a nontrivial dependence 
on t\/t 2 , which does not reduce to a simple power p8[ . 
It follows that the putative scaling law (|l6| ) will fail at 
0(e 2 ): The dynamics of the global order parameter are 
non-Markovian in general, and the exponent 9 is an in- 
dependent critical exponent. Similar conclusions follow 
consideration of the next term in the 1/n expansion. 

The exponent 9 was measured numerically for 2-d Ising 
systems, using a finite-size scaling technique for square 
lattices of linear size L = 24, 32, 48, 64, 96 and 128, with 
periodic boundary conditions. Starting from a random 
initial condition, the systems were evolved using heat- 
bath Monte-Carlo dynamics at the bulk critical coupling 
K c = [m(l + \/2)]/2. Each system was evolved until the 
global magnetization first changed sign. The fraction p(t) 
of surviving systems was then computed at each time t, 
over a number of runs varying from 200 000 for L = 24 
to 91 008 for L = 128. Finite-size scaling suggests the 
scaling form p(t) = t- e f{t/L z ) = L~ 9z f(t/L z ), where 
z is the dynamic exponent. We therefore plot L 0z p(t) 
against t/L z , and vary 9 for the best data collapse. The 
dynamic exponent was taken to be z = 2.15 (HI. Data 
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for t < 20 were discarded. The best collapse was ob- 
tained for 9z = 0.50 ± 0.02. Scaling plots for 6z = 0.48, 
0.50 and 0.52 are shown in Figure 1. We have tried val- 
ues of z other than 2.15, but find that z — 2.1 and 2.2 
give significantly worse data collapse. 

A finite-size scaling analysis is essential here, as the 
data show significant curvature in the 'early' time regime, 
even for the largest systems studied (128 2 ). In the scal- 
ing form L Sz p(t) = F(t/L z ), the scaling function F(x) 
must vary as x~ e for small x, but the 'small- a;' regime in 
the data is not extensive enough to extract the exponent 
from this part of the plot alone. In the large- a; regime, 
one expects F(x) ~ exp(— const x), since L z is the char- 
acteristic relaxation time of the system. This behavior 
is confirmed in studies of smaller systems, where longer 
runs arc feasible. The final part of the scaling plots in 
Figure 1 show the entry into this exponential regime. 



an independent critical exponent, not related by scaling 
laws to other critical exponents, although the relation 
( plf ) is exact for n — oo, to first order in e = 4 — d, and 
for the d — 1 Ising model (for which the dynamics are 
Markovian) , and seems to be numerically quite accurate 
for the d — 2 Ising model. The corresponding exponent 
for the global order parameter following a quench into the 
ordered phase is also of interest, and is currently under 
investigation by numerical simulations. 
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SISSA (Italy), and Amos Maritan in particular, for hos- 
pitality during the final stages of this work, and Alan 
McKane for discussions. 
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FIG. 1. Finite-size scaling plots for the 'persistence' p(t) 
(the fraction of systems whose total magnetization has not 
changed sign) for the d = 2 Ising model at T c , with z = 2.15 
and a = 6z = 0.48, 0.50 and 0.52. For clarity, the 0.50 and 
0.52 data have been moved up by 0.2 and 0.4 respectively. 

It is interesting to compare the numerical result for 
9z with the prediction of the 'scaling law' (16). Using 
A ps 1.59 ± 0.05 [§, and the exact result r) = 1/4, ^) 
gives Oz — 0.47±0.05, consistent with the measured value 
0.50 ±0.02. This suggests that non-Markovian violations 
of the relation ( |l6| ) may be quite small in practice. 

In summary we have identified a new exponent 9 for 
critical dynamics. It is the analog of the 'persistence' ex- 
ponents discussed in a number of other contexts recently, 
and characterizes the time-dependence of the probability 
that the global order parameter has not changed sign up 
to time t after a quench to the critical point from the 
disordered phase. We have argued that 9 is in general 
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